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In a recent paper ([l] = quant-ph/0606035 ) it is shown how the optimal recovery operation in 
an error correction scheme can be considered as a semidefinite program. As a possible future 
improvement it is noted that still better error correction might be obtained by optimizing the 
encoding as well. In this note we present the result of such an improvement, specifically for the 
four-bit correction of an amplitude damping channel considered in [1] . We get a strict improvement 
for almost all values of the damping parameter. The method (and the computer code) is taken from 
our earlier study of such correction schemes (quant-ph/0307138). 
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I. INTRODUCTION 

In a recent paper Fletcher, Shor and Win 1] analyze 
error correction schemes obtained by ab initio optimiza- 
tion, rather than the adaptation of classical coding tech- 
niques. The basic idea H, H, ^ IE 111 is that both en- 
coding £ and recovery (or decoding) TZ can be arbitrary 
channels, and that for a given number n of invocations 
of a noisy discrete memoryless channel the objective is 
to bring the channel f A/"®"??, as close to the identity as 
possible. If the so-called channel fidelity used as a fig- 
ure of merit (as in and 0) the optimization is clearly 
a semidefinite problem, by virtue of Jamiolkowski-Choi 
duality 

In |2| we presented an alternative algorithm for this 
semidefinite problem, and used it to generate optimal 
codes for various noisy channels, by alternatingly opti- 
mizing the encoding channel and the recovery channel. In 
0] only the recovery is optimized, which already gives a 
marked improvement some over previously known codes, 
specifically for the case of the four-bit correction of the 
amplitude damping channel 8] . The possibility of further 
improvements by optimizing also the encoding is noted 
in the discussion (citing also Q)- 

It so turns out that the required computation (even 
for the same test case) was already done in the autumn 
of 2003 in a collaboration between the first two authors 
and the third author of this note, with the aim of check- 
ing the power-iteration method of ^ against the better 
established semidefinite method. Since these results di- 
rectly support the perspective forwarded in , we felt it 
appropriate to make them immediately available. 



II. RESULTS 

For the theoretical background we refer to either or 
0. In both papers similar ideas and notations are used, 
so they should be readily accessible from each other. The 
amplitude damping channel is the qubit channel J\f — M-y 
with Kraus operators 
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FIG. 1: Error correction results for the amplitude damping 
channel (four copies) using three methods: (a) no coding, (b) 
optimized decoding by Fletcher et al. (with encoding by 
Leung et al. Q) and (c) our iterative optimization of both, 
encoding and decoding Q]. 



These channels form a semigroup [N-yNe ~ J^-ye), which 
contracts to the first "unexcited" basis state. The chan- 
nel fidelity of a noisy channel J\f is defined as 



F{J\f) := (f^|(A/'(»id)(|f^)(f^|)|f^), 



(2) 



Kn = 



1 




and K,^(^Q^ 



(1) 



where 51 is a maximally entangled vector. Note that this 
requires input and output of the channel to be systems 
with the same Hilbert space, which is adapted to com- 
paring the channel with the identity, which is the unique 
channel with F{Af ) = 1. This is also closely related [13 
to the average fidelity for pure input states (with the 
average taken according to the unitarily invariant mea- 
sure) . The main virtue of choosing this fidelity as a figure 
characterizing the deviation from the identity is that it 
is linear in A/". Such a linear criterion is possible only 
because the ideal channel is on the boundary of the set 
of channels. 

The curve 7 1-^ F(A/'-y) is the dotted line in Fig. 1. 
The other lines represent 7 1-^ F{£^N!^'^TZ^) for various 
choices of £^ and TZ^ . The dashed line uses the encoding 
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5^ by Leung et al. 

|0)i- -^(10000) + 11111)), 

|l)L = -i=(|0011) + 11100)), 

with an optimized decoding TZ^ 0. In fact, we computed 
this Hne by our routines, and it coincides to within pixel 
resolution with the graph in 0. The solid line is the 
result of the iteration j^], in which £^ and TZ^ are al- 
ternatingly optimized, keeping the other operation fixed. 
This iteration has a strict improvement over the Leung 
code for < 7 < 1. 

III. DISCUSSION 

The methods in and have the following charac- 
teristic features: 

• For a known channel these methods yield excellent 
results, without using any special properties of the 
channel like symmetry etc.. 

• The optimization of either encoding or decoding is 
a semidefinite problem for which the solution is a 
certified global optimum. The process of alternat- 
ingly optimizing these therefore improves the ob- 



jective in every step, and hence converges to a lo- 
cal optimum. However, there is no guarantee for 
having found the global optimum. 

• The methods suffer from the familiar explosion of 
difficulty as the system size is increased. Correc- 
tion schemes like the five qubit code can still be 
handled on a PC, but a nine qubit code would in- 
volve optimization over 2^" x 2^'^-matrices, which 
is set up by multiplying and contracting some ma- 
trices in 2^* dimensions. This may be possible on 
large machines, but it is clear that these methods 
are useless for asymptotic questions. 

• The iteration method replacing the semidefinite 
package in j^] has a slight advantage here, because 
it works with a fixed number of Kraus operators. 
So for the encoding one can put in by hand an iso- 
metric encoding, which, as our study shows, is often 
optimal. This cuts down on dimensions, at least for 
the optimization of encodings. 

• For asymptotic coding theory one still needs codes, 
which can be described also for very large di- 
mensions, be it by explicit parameterization or by 
a characterization of typical instances of random 
codes. It is here that methods transferred from 
classical coding theory will continue to be useful. 
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